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ABSTRACT 

We report on a detailed analysis of the so-called ~1 Hz quasi-periodic oscillation (QPO) in the eclipsing and 
dipping neutron-star low-mass X-ray binary EXO 0748-676. This type of QPO has previously been shown 
to have a geometric origin. Our study focuses on the evolution of the QPO as the source moves through the 
color-color diagram, in which it traces out an atoll-source-like track. The QPO frequency increases from ~0.4 
Hz in the hard state to ~25 Hz as the source approaches the soft state. Combining power spectra based on QPO 
frequency reveals additional features that strongly resemble those seen in non-dipping/eclipsing atoll sources. 

We show that the low-frequency QPOs in atoll sources and the ~1 Hz QPO in EXO 0748-676 follow similar 
relations with respect to the noise components in their power spectra. We conclude that the frequencies of both 
types of QPOs are likely set by (the same) precession of a misaligned inner accretion disk. For high-inclination 
systems, like EXO 0748-676, this results in modulations of the neutron-star emission due to obscuration or 
scattering, while for lower-inclination systems the modulations likely arise from relativistic Doppler boosting 
and light-bending effects. 

Keywords: accretion, accretion disks - stars: neutron - X-rays: binaries - X-rays: individual (EXO 0748-676) 


1. INTRODUCTION 

Low-frequency quasi-periodic oscillations (LF-QPOs) are 
a common phenomenon in the X-ray light curves of low- 
mass X-ray binaries. In power spectra they show up as nar¬ 
row peaks with frequencies ranging from ~1 to ~70 Hz in 
neutron-star low-mass X-ray binaries (NS-LMXBs) and from 
^0.01 to ^30 Hz in black-hole low-mass X-ray binaries (BH- 
LMXBs). Comparisons of the power spectra of NS-LMXBs 
and BH-LMXBs suggest that the most common types of LF- 
QPOs in these systems (as well as some of the other variability 
components) may have a similar origin (Wijnands & van der 


|Klis|1999[|Psaltis et al.|1999lfBelloni et al. 2002) 

For the most common type of LF-QPOs in BH-LMXBs, 
the so-called type-C QPOs ( Casella et al.||2005) , there is in¬ 
creasing evidence that points towarcfa geometric origin: these 
QPOs are stronger in sources that are viewed at higher in¬ 
clinations ( |Schnittman et al.||2006| |Heil et aT7||20 1 5[ |Motta| 
et al.|2015 i, and the strength of the Fea emission line in GRS 


T9f‘5+TU5~has been found to be strongly modulated as pa rt of 
the type-C QP O cycle ( |Miller & Homan|2005| |Ingram & van| 
|der K lis 2015). Both these results are consistent with the QPO 
being generated b y a m isaligned and precessing inner accre¬ 
tion disk (Schnittma n et al. |2006| l, with relativistic Lense- 
Thirring precession being a commonly suggested mechanism 

(e.g., [Ingram et al.|2009| >. 

In seven dipping and/or eclipsing NS-LMXBs a type of LF- 
QPO has been observed whose properties suggest (even more 
directly) that it also has a geometric origin (Homanj |2012| i. 
These QPOs are often referred to as the “~1 Hz QPOs”7 due 
to the fact that they were initially discovered around 1 Hz in 
three different sources ( [Homan et aL|1999}|Jonker et al.| 1999| 
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20()()i). The geometric nature of these QPOs was deduced 
from the fact that these QPOs persist during type I (thermonu¬ 
clear) X-ray bursts with a fractional amplitude that is consis¬ 
tent with that observed during the persistent (non-burst) emis¬ 
sion, and that their strength is nearly independent of energy. 
Both these properties suggests that (part of) the emission from 
the neutron star is obscured or scattered out of the line of sight 
by the accretion flow. Homan (2012) recently suggested that 
an inclined and precesling inner accretion flow might be a 
possible explanation for the observed properties of the ~1 Hz 
QPOs. The misalignment angle between the inner flow and 
binary plane would have to be ~15°-20° for this type of QPO 
to be observed only in dipping and/or eclipsing NS-LMXBs. 
The relation between these ~1 Hz QPOs and the LF-QPOs 
observed in the non-di pping/ec lipsing NS-LMXB is not clear 
at this point. Although Homan]( [201 2\ pointed out some simi¬ 
larities (e.g., both being mainly observed in the harder spectral 
states) the ~1 Hz QPOs in the dipping/eclipsing sources have 
not been studied in enough detail to make a proper compari¬ 
son. 

In this paper we investigate the ~1 Hz QPO observed in 
the dipping and eclipsing NS-LMXB EXO 0748-676. Of the 
seven sources in which ~1 Hz QPOs have been observed, 
EXO 0748-676 has the largest reported frequency range. It 
also has the highest number of Rossi X-ray Timing Explorer 
(RXTE ) observations and is known to make regular state tran¬ 
sitions. We focus on how the ~1 Hz QPO evolves along 
the tracks traced out in the color-color diagram (CD; Section 
ED and perform a comparison of the observed power spec¬ 
tra with those of non-dipping/eclipsing NS-LMXBs (Section 
|3.2| ). Our results are discussed in Section[4] 

2. OBSERVATIONS AND DATA REDUCTION 


The observations of EXO 0748-676 analyzed in thi s paper 
were made with the Proportional Counter Array (PCA; [Jahod a 
|et al.|[2006) > onboard RXTE ( |Bradt et al.||1993| l. EXO 0748- 
676 was observed 752 times with the PCA during the lifetime 
of RXTE. A number of observations had to be discarded for 
various reasons: they were made after the source had returned 
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Figure 1. Long-term light curves of EXO 0748-676. In the top panel we 
show 7-day averages from the RXTEI ASM 1.5—12 keV). In the bottom 
panel we show data from the RXTEfPCA (~2-60 keV), with each data point 
representing the average count rate of a single group. Errors on the PC A 
count rates are smaller than the size of the plot symbols. Around MJD 54700 
(2008 August) the source returned to quiescence. 


to quiescence ( jWolff et al.||2008b|a| ), important data modes 
were missing, or not enough data were left to produce power 
spectra of the required length (128 s) after the recommended 
filtering criteria had been applied. This left 718 useful ob¬ 
servations. To increase the sensitivity to QPOs, observations 
were combined into groups. Consecutive observations were 
added together if the time gap between them was less than 
12 hr, resulting in 203 groups of observations. In all almost 
all cases these groups covered periods less than a day, but for 
seven groups the period was longer (up to ~2 days). These 
seven groups were split into two or three shorter segments, 
covering less than a day. The final number of groups was 211, 
with individual groups containing between 1 and 11 observa¬ 
tions. The group exposure times ranged from ~350 s to ~40 
ks, with an average of ~8 ks. Note that this process of com¬ 
bining observations is justified by the fact that the ~1 Hz QPO 
in EXO 0748-676 was seen to vary in frequency by less than 
~15% on a timescale of a day (Homan et al .{1999) . 

We used standard 2 mode data to create light curves, 
hardness-intensity diagrams (HIDs), and CDs. Data from all 
available Proportional Counter Units (P CUs) were used. Fol¬ 
lowing the procedures laid out in Fridriksson et al.| (f2015 ), 
count rates from each PCU were corrected for changes rn 
the response, using long-term light curves of the Crab, and 
then normalized to those of PCU 2. We define a soft color 
as the net counts in the 4.0-7.3 keV band divided by those in 
the 2.4-4.0 keV band, and a hard color (hardness) as the net 
counts in the 9.8-18.2 keV band divided by those in the 7.3- 
9.8 keV band. The intensity we use for the HIDs is the net 
count rate per PCU in the 2-60 keV band. The standard 
2 data were barycenter corrected and eclipses were removed 
using the ephemeris from |Wolff et al.|(|2009). Type I X-ray 
bursts were removed using standard 1 mode data, fol¬ 


lowing the procedures laid out in Remillard et al. 

and high-time-resolution data were fil- 


The standard 2 

tered following a similar procedure. 
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Figure 2. Hardness-intensity diagram (left) and color-color diagram (right) 
of EXO 0748-676. Each data point represents the average of a single group. 
See Section[2]for definition of colors. 

removed from the standard 2 data, since they strongly af¬ 
fect the average count rates and colors of the observations. 
The removal of dips had to be done manually and was based 
on visual inspection of the light curves and soft-color curves. 
Given the strong differences between observations in terms 
of dipping activity it was difficult to remove dips in a fully 
consistent manner for the entire data set. We did not remove 
the dips from the high-time-resolution data used for the cre¬ 
ation of the power spectra, since the ~1 Hz QPO can still be 
detected during dips (Homan et al.j| 1999j l. 

Power spectra were created from light curves extracted 
from the high-time-resolution data modes. Light curves were 
created with a time resolution of 1/8192 s (except for the first 
observation, 10068-03-01-00, which had a maximum time 
resolution of 1/2048 s), covering absolute PCA channels 0- 
35. This channel band effectively corresponds to an energy 
range of ~ 1.5-9.5 keV during Gain Epoch 1, increasing to 
~2.1—15 keV during Gain Epochs 4 and 5. We note that more 
than 600 of our observations were made during Gain Epochs 4 
and 5. To create the power spectra we performed Fast Fourier 
Transformations of 128 s light-curve segments. The power 
spectra were averaged for each group and rms normalized 
(Belloni & Hasinger||1990| |Miyamoto et al.||1991| |van der| 
k.lis||1995| l. These normalized power spectra were rebinned 
logarithmically to 70 bins per frequency decade; for plotting 
purposes the power spectra were sometimes rebinned even 
further. Prior to fitting we subtracted the deadtime-modified 
Poisson level, which we approximated by taking the aver¬ 
age power in an interval of a few 100 Hz centered around 
~2()00 Hz. While this method may introduce some uncertain¬ 
ties in the properties of the high-frequency features, it should 
have a negligible effect on the low-frequency features in the 
power spectra (which are the focus of this paper). The result¬ 
ing power spectra were fitted up to 2000 Hz with a combina¬ 
tion of a power law (P(v) oc ^ - “), and one or more Lorentzians 
(P(v) = (r 2 A/7r)[A 2 + (^-^o) 2 ] -1 )- Here vq is the centroid 
frequency, A the half-width-at-half-maximum, and r the in¬ 
tegrated fractional rms (from -oo to oo). Instead of vq and 
A we will quote the frequency z/ rnax at which the Lorentzian 
attains its maximum in vP{v) and the quality factor, Q, where 
zzq( 1 + 1 /4Q 2 ) 1,/2 and Q = zzq/2A (feelloni et al. 2002). 


The fractional rms amplitudes quoted in tins paper were the 
integrated power between 0 and oo Hz for the Lorentzians, 
and between 0.01 and 100 Hz for the power-law component. 
Errors on fit parameters were determined using A\ 2 = 1 (i.e. 
68% confidence). 


Absorption dips were 
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Figure 3. Color-color diagram of EXO 0748-676, showing groups with ~1 
Hz QPO detections (red) and those without (gray). The black line shows the 
spline used to calculate the S a values for each data point. The two large black 
circles indicate the locations of S a = 1 and S a = 2. The smaller black circles 
represent steps in S a of 0.1. Error bars are omitted for clarity. 


3. ANALYSIS & RESULTS 
3.1. The ~1 Hz QPO in EXO 0748-676 

In Figure [I] we show long-term light curves of EXO 0748- 
676 from the RXTE All-Sky Monitoif] (ASM; Levine et al. 
1996) and the PC A. The ASM light curve shows 7-day 


erages, while the PCA light curve shows the average count 
rate for each group. Especially during the second half of 
RXTE’s lifetime EXO 0748-676 was quite variable, with oc¬ 
casional several-week-long outbursts. In Figure [2] we show 
the HID and CD corresponding to the PCA light curve. The 
patterns traced out in those diagrams are very similar to those 
of some of the low-l uminosit y (atoll) N S-LMXBs (see, e.g., 
[Gladstone et al.| ( |2007| ) and |Fridriksson| ( |2011| ) for colle ctions 
of at oll-source CDs/HI Ds). As a lready pointed out by |Ponti| 
|et al.||2014| ) and |Munoz-Darias et al. ( |2Q 14] ), hard and soft 
states can be identified in the CD/HID of EXO 0748-676. In 
the HID in Figure [2] we label the hard, intermediate, and soft 
states; these states are often also referred to as the extreme 
island state, the island state and the (lower) banana state, re¬ 
spectively. A comparison between the CD of EXO 0748-676 
and those of the atoll sources presented in |Fridriksson| (j2011 ) 
indicates that the soft color of EXO 0748-676 spans an un¬ 
usually large range. This is likely the result of incomplete dip 
removal, which mostly affects that low energy bands. We note 
that before the dip removal the soft color range was ~ 1.3-5.1. 

Based on a visual inspection of all 211 power spectra for 
the presence of a ~1 Hz QPO, we selected 171 for further 
analysis. Following the analysis of Homan et ah (1999) these 
171 power spectra were initially fit with a combination of 
a power law and a single Lorentzian. While in most cases 
this model resulted in acceptable fits (average xj? ed ~ 1.18), in 
some power spectra the QPO profile was not accurately fit be¬ 
cause of a broad underlying noise component. Adding a zero- 
centered Lorentzian to the fit function solved this issue and 
resulted in an average xjr ed of ~ 1.05. QPOs with single-trial 



Figure 4. Properties of the ~1 Hz QPO in EXO 0748-676 as a function 
of S a . From top to bottom we show frequency, Q- value, and fractional rms. 
The errors on the frequency are typically smaller than the symbol size. The 20 
gray horizontal bars in the upper panel show the boundaries between intervals 
used for grouping power spectra based on QPO frequency (see main text). 


significance*^] larger than 3cr were detected in the power spec¬ 
tra of 157 groups, with frequencies ranging from ^0.43 to 
~25.5 Hz. Together these 157 groups represent 558 (~78%) 
of the 718 observations we analyzed. 

To see in which states the ~1 Hz QPO was detected we 
show a color-coded CD in Figure [3} the groups with LF-QPO 
detections are colored red, and those without are colored gray. 
It is clear from this figure that the ~1 Hz QPO was predomi¬ 
nantly detected in the hard state, with a few detections in the 
intermediate state. In order to study how the QPO properties 
evolve along the track in the CD, we parameterized the posi¬ 
tion along the track by manually drawing a spline onto which 


all data points were projected (see, e.g.. Dieters & van der 
Klis 2000). To correct for the differences in scale between the 
horizontal and vertical axes in Figure [3] the soft-color values 
were scaled down by a factor of 0.286 during the projection 
procedure, thus ensuring that the data points were projected 
onto the nearest spline point in the figure. Two normal points, 
depicted by the large black circles in Figure[3j were chosen to 
set the scale for the track parameter S a . .S' a = 1 was chosen to 
fall close to the end of the dense cluster of points in the hard 
state of EXO 0748-676, while S a = 2 was chosen to coincide 
with the boundary of the intermediate state and the soft state 
(most clearly seen in the HID in Figure [2] hard color ~0.8). 
The smaller black circles represent steps of 0.1 in S a . Using 
this parametrization we find that no ~ 1 Hz QPO was detected 
significantly for S d values greater than 1.73. Note that the ob¬ 


servation th at showed the 695 Hz QPO reported by Homan & 
| van der Klis| < ]2000| ) has an ,S' a value of ~2.1. 


4 Public ASM data were obtained from http: / /xte . mit. edu 


5 The integrated power divided by its negative error. 
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Figure 5. Three representative power spectra of EXO 0748-676 from 
frequency-selected groups. Fits to the power spectra are shown as solid 
black curves, while individual components are shown in gray. The various 
power-spectral components are identified in the top panel. The ~1 Hz QPO 
is marked as Llf- 

In Figure [4] we show the QPO frequency, Q-value and rms 
amplitude as a function of S. d . The QPO frequency shows a 
clear increase with .S' a , although there is considerable scatter 
at low .S'a values. The (Lvalue shows no clear dependence on 
,S a , while the rms appears to increase toward higher ,S a . 

Part of the scatter observed in Figure |4]is likely caused by 
the fact that we were not able to perfectly remove all dipping 
intervals from the data. Since dipping is most prevalent in the 
hard state (i.e., at low .S a ) and increases the soft color, the re¬ 
moval of dips will move observations parallel to the spline in 
the CD, towards higher .S' a . Various gradations of incomplete 
dip removal will therefore result in a range of S a values for ob¬ 
servations whose non-dipping emission would have the same 
.S’a value. To reduce the impact from these effects, and also to 
boost the signal-to-noise ratio of our power spectra, we per¬ 
formed additional grouping of our observations. Observations 
from groups with similar QPO frequencies were combined 
into 21 new groups. At the low-frequency end the width of the 
frequency selections for the new groups was 0.1 Hz, increas¬ 
ing to larger values for higher frequencies; the gray horizon¬ 
tal lines in the top panel of Figure [4] indicate the boundaries 
separating the 21 frequency-based groups. These new groups 
contained between 1 and 17 of the old groups, and between 
1 and 69 of the original observations. To determine ,S' a values 
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Figure 6. Properties of the ~1 Hz QPO in EXO 0748-676 as a function of 
S a . From top to bottom we show frequency, Q- value, and fractional rms. The 
errors on the frequency are typically smaller than the symbol size. 

for each of the new groups we reprojected the averaged CD 
data points onto the previously used spline. 


As a result of the increased signal-to-noise we were able to 
detect additional features in the power spectra, especially for 
the power spectra at low .S' a . In Figure [5] we show three rep¬ 
resentative power spectra, in which the — 1 Hz QPO has fre¬ 
quencies of —0.5 Hz, —1.4 Hz, and —6.6 Hz (top, middle, and 
bottom panel, respectively). Fits to these power spectra are 
shown as well, including the individual components (a power 
law and up to five Lorentzians). As can be seen, in all cases 
the —1 Hz QPO is accompanied by several broad components 
and a power-law noise component. Base d on a comparison 
with other NS-LMXBs (see Section |3.2| i and following the 
nomenclature used in prior works, we identify the following 
components in the power spectrum in the top panel of Fig¬ 
ure [5] a break (Lb), the —1 Hz QPO, which we identify as a 
LF-QPO (Llf), a ‘hump’ (LQ, and two broad features that are 
commonly associated with the lower and upper kHz QPOs (L\ 
and L u ). Given the quality of the power spectra, the lower and 
upper kHz QPOs are typically not well constrained and they 
will not be discussed in detail. The fit parameters of the —1 
Hz QPO (Llf) and the two low-frequency noise components 
(Lb and Lh) are given in Table [I] The power-law component 
gradually steepened from a wl.5 at low QPO frequencies 
to a k 2.5 at high QPO frequencies, while it decreased in 
strength from —7%-8% to —2%-3%. 

In Figure [6] we once again show the evolution of the — 1 
Hz QPO properties as a function of S a , this time using the 
measurements from the frequency-selected groups. As can 
be seen from the top panel, combining groups based on fre- 
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Table 1 

Power Spectral Fit Parameters 





Llf 



L b 



L h 


No." 

Sa 

^max 

(Hz) 

Q 

rms 

(%) 

k'max 

(Hz) 

Q 

rms 

(%) 

^max 

(Hz) 

Q 

rms 

(%) 

1 

1.142±0.003 

0.483±0.003 

2.9±0.2 

9.8±0.3 

0.20±0.11 

0.83±0.19 

6.5±0.7 

1,75±0.13 

0.58±0.19 

9.8±1.3 

2 

1.116±0.003 

0.574±0.004 

3.1 ±0.3 

8.9±0.4 

0.19±0.16 

0.7±0.3 

5.9±0.9 

2.1±0.3 

0.2±0.2 

12.9fjj 

3 

1.122±0.003 

0.665±0.003 

3.5±0.2 

11.1 ±0.3 

0.27±0.19 

0.8±0.2 

6.7±0.8 

2.11±0.16 

0.34±0.15 

14.8± 1.3 

4 

1.158±0.001 

0.749±0.003 

3.14±0.17 

10.9±0.3 

0.27±0.16 

0.8±0.2 

5.1±0.6 

2.25±0.10 

0.42±0.09 

15.8±1.1 

5 

1.090±0.002 

0.857±0.003 

3.36±0.16 

10.2±0.2 

0.31±0.15 

1.7±0.4 

3.5±0.4 

2.66±0.18 

0.23±0.07 

16.8± 1.0 

6 

1.143±0.001 

0.952±0.004 

3.4±0.2 

12.3±0.4 

0.3±0.2 

1 4+0.7 

1 ’ -0.4 

4.0±0.7 

2.8±0.2 

0.30±0.08 

18.2± 1.2 

7 

1.146±0.002 

1.040±0.004 

3.5±0.19 

10.2±0.3 

0.40±0.03 

1 4 + 0.7 

1 ’ -0.4 

3.3±0.7 

3.2±0.3 

0.23±0.10 

15.7± 1.7 

8 

1.200±0.002 

1.182±0.007 

3.1 ±0.2 

10.2±0.4 

0.50±0.05 

0 9 +0 - 9 
U - V -0.5 

4.3± 1.3 

3.4±0.3 

0.38±0.19 

13.9+™ 

9 

1.175±0.002 

1.397±0.007 

3.25±0.19 

10.7±0.3 




4.5±0.5 

0.04±0.08 

17.6±0.9 

10 

1.167±0.002 

1.598±0.009 

2.92±0.18 

11.6±0.3 




5.7±0.7 

0* 

18.3±0.8 

11 

1.209±0.002 

1.901±0.013 

3.1 ±0.3 

10.3±0.4 




4.5±1.2 

0.09±0.12 

1 5 3 +1 ' 6 
aj - :) -0.5 

12 

1.261±0.002 

2.430±0.014 

2.76±0.19 

13.2±0.5 

1 09+ 047 

0.6±0.3 

5.7±1.9 

7.4±0.7 

0.3±0.2 

15.4± 1.4 

13 

1.298±0.005 

2.84±0.05 

2.1 ±0.3 

11.0±0.7 




7 1 +L8 
'■ -1.2 

0* 

13.3± 1.0 

14 

1.303±0.003 

3.75±0.03 

2.9±0.3 

11.6±0.5 




8.0±1.0 

0 b 

14.9±0.9 

15 

1.322±0.005 

4.61±0.05 

2.9±0.3 

12.5±0.6 




7.6±1.2 

0* 

15.3± 1.0 

16 

1.361±0.009 

5.43±0.11 

2.9±0.6 

12.7± 1.1 




15+f 

0* 

13 1 +2 - 5 

17 

1.370±0.003 

6.68±0.10 

1.8±0.3 

14.7±0.9 

2.8±0.9 

1 o +0 - 8 

1 ’-0.4 

5 6 +2 - 5 
* 3 ’ —1.3 

20±3 

i o +0 - 9 

A ■ - 0.6 

9 1 +2 - 9 
.3 

18 

1.407±0.004 

8.96±0.16 

2.3±0.4 

12.9±0.9 




16±3 

0* 

18.0±1.1 

19 

1.516±0.005 

12.8±0.2 

2.5±0.5 

12.0±1.0 




22±8 

0* 

13 l +2 -° 
AJ-A -0.9 

20 

1.631 ±0.009 

17.0±0.4 

4.6+i 6 5 

9.4±1.5 




28^2 

o b 

!3.6^;| 

21 

1.729±0.015 

25.3L1.2 

3.5±1.6 

11.3± 1.7 








a Number of frequency-selected group, 
b Value was fixed. 


quency greatly reduces the scatter at low S a compared to Fig- 
ure[4] although some scatter is still present. Also, for the same 
reason no data points are present any longer below S d = 1.09. 
As before, the QPO frequency shows a strong increase with 
S d . The (Lvalue decreases slowly up to S a « 1.4-1.6 and ap¬ 
pears to increase after that. Although the fractional rms am¬ 
plitude of the QPO still shows considerable scatter at low S d 
values, overall there is an increasing trend up to S d ~ 1.4. The 
rms appears to drop again toward higher S d values. 

The Lh noise component is detected in all frequency- 
selected power spectra, except for the last one. The Lb noise 
component is consistently detected up to QPO frequencies 
of ~1.2 Hz. At higher QPO frequencies the L h component 
becomes broader and as a result the Lb component becomes 
more difficult to detect. It is only detected twice more, when 
the Lh component is narrower. 

3.2. Comparison with other Atoll Sources 

Since EXO 0748-676 shows atoll-like patterns in its CD 
and HID, we investigated whether the ~1 Hz QPO in EXO 
0748-676 is possibly related to the LF-QPOs seen in the 
non-dipping/eclipsing atoll sources. Due to their relatively 
low rms amplitudes (~few percent) detections of LF-QPOs in 
atoll sources are not very common (especially compared to the 
LF-QPOs in the more luminous Z sources). We first searched 
the literature for RXTE observations of non-dipping/eclipsing 
atoll sources with LF-QPOs that have frequency ranges (par¬ 
tially) overlapping that of the ~1 Hz QPO in EXO 0748-676. 
Several examples were found: 4U 1636-53 (Altamirano et al. 
|2008a] >, 4U 1608-52 dvan Straatenetal.|2003| ), IE 1724-3045 


(Altamirano et a l.[2008b| , and 4U 1820-30 (Altamirano et ah 


2005|l. In addition we also analyzed observations of 


U1812— 

12 and found LF-QPOs around 0.3 Hz. In Figure [7 
pare the power spectra of the first three of these atoll sources 
with three power spectra of EXO 0748-676. 

For similar QPO frequencies we find that the power spec¬ 
tra of EXO 0748-676 and the three atoll sources are to a 
large extent similarly shaped, with the accompanying noise 
and (broad) QPO components having similar frequencies as 
well. Most importantly, like the ~1 Hz QPO in EXO 0748- 
676, the LF-QPO component in the three atoll sources has a 
frequency between that of the Lb and Lh noise components (al¬ 
though in IE 1724—3045 it is closer to the Lh component than 
in the other sources). However, there are several differences 
as well. The ~1 Hz QPO in EXO 0748-676 is considerably 
stronger than the LF-QPOs detected in the atoll sources (in 
the same frequency range). As can be seen in Table[l] the ~ I 
Hz QPO in EXO 0748-676 has a fractional rms amplitude 
in the range ~9%-15%. The atoll source LF-QPOs, on the 
other hand, have fractional rms amplitudes of ~2%-5%. On 
average the (Lvalues of the atoll sources LF-QPOs tend to be 
higher than those of the ~1 Hz QPO in EXO 0748-676; this is 
especially clear for the examples shown in Figure IT] In EXO 
0748-676 the (Lvalues ranged from ~1.8 to ~4.67 while in 
the atoll source values of ~ 1.6—8.8 are found (see above ref¬ 
erences). While grouping of observations may have led to 
some broadening of the QPO in EXO 0748-676, similar Q- 
values were already reported for QPO detections in individual 
observations by Homan et al. (jl999 1 . 
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Frequency (Hz) Frequency (Hz) Frequency (Hz) 

Figure 7. A comparison of power spectra of EXO 0748-676 with those of three atoll sources. The following RXTE ObsIDs were used for the atoll sources: 
60032-05-21-00(0) for 4U 1636-53, 30062-02-02-00(0) for 4U 1608-52, and 10090-01-01-000 for IE 1724-3045. 


There are some differences in the behavior of the other 
power spectral components as well. The power spectra of 
EXO 0748-676 contain a clear power-law noise component, 
which is absent in the atoll source power spectra. Although 
some atoll sources do show such a power-law component 
in their power spectra, it typically does not appear until the 


source is closer to the soft state (see, e.g., van Straaten et al. 
2003}. We note that strong power-law components are present 


in some of the other dipping and/or eclipsing sources as well 
(Jonker et al. ]1999[[2000| >, suggesting its presence is related to 
viewing angle. Another difference pertains to the two broad 
high-frequency components (which we tentatively identified 
as the lower and upper kHz QPOs); these tend to be stronger 
in EXO 0748-676 than in the other atoll sources, as can be 
seen from Figure IT] 

As the ~1 Hz QPO in EXO 0748-676 and LF-QPOs in 
the atoll sources appear to have similar frequencies relative 
to those of the two main (Lb and L h ) noise components, it 
is interesting to explore whether the QPOs also follow simi¬ 
lar relations in QPO frequency vs. noise frequency diagrams. 
Such studies have been performed for various power-spectral 


components in NS-LMXBs and BH-LMXBs (see, e.g., Psaltis 

et al. 1999 

Wijnands & van der Klis 1999; 

Belloni et al. 

2002). In Figure 18 

EXO 0748-676 and 

we plot the frequency of the QPOs in 
several atoll sources as a function of the 


frequency of the break component Lb deft panel) and as a 
function of the frequency of the hump component Lh (right 
panel). In addition to the atoll sources mentioned earlier, we 
also include data from the dipping NS-LMXB EXO 1745- 
248 (stars), which also shows a ~1 Hz QPO ( |Homan|2012| >. 
For this source the break frequency was poorly constrained. 
Finally, for reference we show in the left panel of Figure[8]the 
original data of Wijnands & van der Klis ( 1999 ) for the Z and 
atoll sources. We noteTTiowever, that for the power spectra of 
the atoll sources Wijnands & van der Klis ( |1999| > mostly used 
the frequency or Lh for that of the QPO. As a result the QPO 
values that we plot for EXO 0748-676 and the atoll sources 
fall well below the relation found by Wijnands & van der Klis 
( 1999] ). 

As can be seen from both panels in Figure [8] EXO 0748- 
676 traces out relations that are largely consistent with those 
of the other atoll sources. We note that IE 1724—3045 (gray 


squares) presents somewhat of an outlier in both panels. The 
QPO frequency for this source lies a factor of ~2 above that of 
the relation traced out by the other sources, suggesting that the 
QPO in IE 1724-3045 might be the second harmonic, with 
the fundamental not being present. In 4U 1812-12 (triangles) 
and EXO 1745-248 (stars) we detect the fundamental as well 
as the second harmonic. Both are shown in Figure [8] with 
the fundamentals falling on the main relation and the second 
harmonics coinciding with the cluster of IE 1724-3045 data 
points. This appears to confirm the harmonic nature of the IE 
1724-3045 QPO. For this reason we also show the QPO fre¬ 
quencies of IE 1724-3045 divided by a factor of 2 in Figure 
[8](red squares); as can be seen, in this way the IE 1724-3045 
data line up well with the main relation. 

4. DISCUSSION 

We have presented a study of the ~1 Hz QPO in the dipping 
and eclipsing NS-LMXB EXO 0748-676. The main goal of 
this work was to understand how this QPO evolves between 
the various spectral states and how its properties compare to 
those of the LF-QPOs in other NS-LMXBs. 

In its HID and CD EXO 0748-676 traces out tracks that 
strongly resemble those of the atoll source NS-LMXBs, with 
hard (extreme island), transitional (island), and soft (banana) 
states. We found that the ~1 Hz QPO in EXO 0748-676 is 
only observed in the hard and intermediate states and that its 
frequency strongly increases as the spectrum softens. After 
combining power spectra based on QPO frequency, a rela¬ 
tively tight relation with the track parameter .S’ a is observed 
(see Figure |6j, with the frequency increasing from ~0.5 Hz 
at Sa ~ 1.2 to ^25 Hz at S 2 « 1.8. This increase in frequency 
from the hard state into the intermediate state, as well as the 
frequency range covered, is similar to what is observed for the 
LF-QPOs in other atoll sources in these sta tes ( [van Stra aten 
|et al.|200'3| | Altamirano et al.|2 005 2008a b). 

The grouping of power spectra from individual ObsIDs 
(first based on time and later also based on QPO frequency) 
resulted in high-quality power spectra that revealed several 
peaked noise components in addition to the QPO and power- 
law noise component already reported in the ~ 1 Hz QPO dis¬ 
covery papers. A comparison of the QPO and noise com¬ 
ponents in the power spectra of EXO 0748-676 with those 

















































































Geometric Nature of Low-Frequency QPOs 


7 


o 

o 


N 

X 


5^ 

O 

c 

(D 

D 

cr 

0 ) 


O 

CL 

a 



N 

X 


>x 

o 

c 

0) 

u 

cr 

0) 


O 

X 

a 


Figure 8. QPO frequency vs. break frequen cy (left p anel) and vs. hump frequency (right panel). See left panel for the meaning of the various symbols. For 
reference we also show the original data from Wijnands & van der Klis 1 1999 1 for ato ll and Z sources. For IE 1724-3045 (gray open squa res) we also sh ow the 
QPO frequencies divided by a factor of 2 (red open squares). Data were taken from|Altamirano et al.||2008al for 4U 1636-53, Altamirano et al. (2005) for 4U 
1820-30, and Altamirano et al.|(200Sb) for IE 1724-3045. For 4U 1608-52 we analyzed KXLE observations 30062-02-02-00(0), tor4U 1812-12 observations 
30701-13-01-|00:02J, and for I NO 1745-248 observation 50054-06-04-03. 


in atoll source power spectra revealed strong overall similari¬ 
ties. Most importantly, like the LF-QPO in atoll sources, the 
~1 Hz QPO in EXO 0748-676 (Llf) is located between the 
break (Lb) and hump (Lh) noise components (Figure [7]) and its 
frequency evolves in a similar way with respect to that of the 
two noise components (Figure[8|. 

These results strongly suggest that the frequencies of the 
~1 Hz QPO in EXO 0748-676 and the LF-QPOs in (other) 
atoll sources are set by the same mechanism. From past work 
( |Homan et al.| 1999| Jonker et al.| [l999[ 2000]) it is clear that 
~1 Hz QPOs in EXO 0748-676 and other dipping/eclipsing 
NS-LMXBs have a geometric nature, in the sense that parts of 
the accretion flow continuously move (rotate) in and out of our 
line of sight to the neutron star, thus modulating the emission 
we see from the neutron star surface and/or boundary layer. 
Homan ( |2012| suggested that a misaligned precessing inner 
flow could provide such a continuously changing geometry. 

Since the non-dipping/eclipsing atoll sources have lower in¬ 
clination angles, and there is no evidence for line of sight 
modulations, such a precessing inner disk would have to 
modulate the X-ray emission in a different manner in those 
sources. Modulation of the X-rays through relativistic beam¬ 
ing and light-bending effects , as suggested for LF-QPOs in 
BH-LMXBs (jSchnittman et al.|2006|), could offer an explana¬ 
tion for the non-dipping/eclipsing sources. We note that these 
effects are likely also at work in dipping/eclipsing sources. 
However, Figure 3 in Schnittman et al. ( |2006 1 suggests that 
the resulting modulations are out of phase with the modula¬ 
tions caused by the line of sight interceptions. This may in 
fact reduce the overall strength of the ~1 Hz QPO, although 
it is still observed to be significantly stronger than the LF- 
QPOs in non-dipping/eclipsing atoll sources. The differences 
in modulation mechanism between the two types of QPOs 
may also explain the (on average) higher rms amplitudes and 
lower quality factors of the ~1 Hz QPOs in comparison to the 
LF-QPOs in the non-dipping/eclipsing atoll sources, although 
detailed simulations are necessary to verify this. In addition to 
resulting in different modulation mechanisms for the QPOs, 
we speculate that the higher inclination of EXO 0748-676 
may also be the cause of the more prominent power-law noise 
component as well as differences in the relative strengths of 
the broad peaked noise components, although the mechanisms 


involved are unclear. 

As mentioned in Section [I] there is increasing evidence for 
a geometric origin of the (type-C) LF-QPOs in BH-LMXBs 
dHeil et al,||2015[ |lngram & van der Kl is |2015i |Motta et al.j 

2015|l, with Lense-Thirring precession (Stella & Vietri|1998 1 


currently being the most promising mechanism (Ingram et al. 
2009|). While we argued above that the ~ 1 H z QPO in EXO 
0748-676 and the LF-QPOs in other atoll sources have a geo¬ 
metric nature as well, it is not clear to what extent their fre¬ 
quencies are set, or dominated, by Lense-Thirring preces¬ 
sion. Classical precession (due to neutron-star oblateness) and 
magnetic precession (due to magnetic torques) may signifi- 
cantly lower the precession frequency (Lai 1999; Shira kawa| 
|& L ai 2002), while radiation feedback from the neutr on star 
could significantly increase the precession frequency ( |Miller| 
|1999| ). The latter mechanism may explain why, for the same 
break (Lb) frequency, the LF-QPOs in the more luminous Z 
sources (the horizontal-branch oscillations) have higher fre¬ 
quencies (by a factor of ~5, see Figure [8j than the ~1 Hz. 
QPO an d LF-QPOs in EXO 0748-6 76 andother atoll sources. 
In fact, [Altamirano et ah] ( |2012] l already pointed out that 
the frequencies of the horizontal-branch oscillations in the 
Z source and 11 Hz pulsar IGR J17480-2446 were too high 
to be explained by pure Lense-Thirring precession, suggest¬ 
ing that radiation feedback can indeed significantly alter (or 
even dominate) the precession frequency in neutron star sys¬ 
tems. However, we stress that the conclusions we have drawn 
about the geometric nature of the ~1 Hz QPO and LF-QPOs 
in atoll sources do not necessarily apply to the horizontal- 
branch oscillations in the Z sources, since these conclusions 
were based on the similarities between the power spectra of 
EXO 0748-676 and the non-dipping/eclipsing atoll sources. 
The horizontal-branch oscillations follow a very different re¬ 
lation in the left panel of Figure [8] (although with a similar 
slope) and it is therefore not clear whether they are the result 
of a similar mechanism as the LF-QPOs in atoll sources. 

In none of the seven NS-LMXBs that have shown ~1 Hz 
QPOs are the QPOs detected in the soft state (Homan 
|2012| . Combined with the fact that LF-QPOs in non¬ 
dipping/eclipsing atoll sources are also rare (or absent) in 
the soft state, this suggests that the precessing part of the 
accretion flow in atoll NS-LMXBs becomes significantly 
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smaller toward the soft state and/or that the misalignment 
angle between the inner flow and neutron-star spin axis de¬ 
creases. Both changes would reduce the effects of line of 
sight interceptions in the dipping/eclipsing systems and of 
the relativistic beaming and light-bending effects in non¬ 
dipping/eclipsing sources. In the case of the ~1 Hz QPO 
a reduction in the geometrical thickness and/or the optical 
thickness of the inner flow may also affect the efficiency of re¬ 
moving photons from the neutron star out of our line of sight 
and therefore reduce the strength of the QPO. Finally, the fre¬ 
quency increase of the ~1 Hz QPOs in EXO 0748-676 and 
the LF-QPOs in the atoll sources indicates that the transition 
radius between the (standard) accretion disk and the precess- 
ing inner accretion flow decreases as the spectrum softens, 
similar to what has been concluded from the type-C QPOs in 
BH-LMXBs ( jlngram & Done[2011| . 

5. CONCLUSIONS 

We have studied the ~1 Hz QPO in the dipping and eclips¬ 
ing NS-LMXB EXO 0748-676 and compared it to the LF- 
QPOs in non-dipping/eclipsing atoll NS-LMXBs. We find 
strong similarities between the two types of QPOs and con¬ 
clude that their frequencies are likely set by the same mecha¬ 
nism. A misaligned precessing inner accretion flow is sug¬ 
gested as the common frequency mechanism; the resulting 
modulations are geometric in nature, but the way in which the 
modulations occur depends on viewing angle. In the (low- 
inclination) non-dipping/eclipsing sources the emission from 
the precessing flow could be modulated through relativis¬ 
tic beaming and light-bending effects, while in the (higher- 
inclination) dipping/eclipsing systems we get an additional 
(and dominating) modulation of the X-rays from the neutron 
star through line of sight interceptions by the precessing inner 
flow. 

J.H. acknowledges financial support from NASA grant 
NNX12AE14G, provided through the Astrophysics Data 
Analysis Program. This research has made use of data ob¬ 
tained from the High Energy Astrophysics Science Archive 
Research Center (HEASARC), provided by NASA’s Goddard 
Space Flight Center. 


REFERENCES 

Altamirano. D., Ingram, A., van der Klis, M., et al. 2012, ApIL, 759, L20 

Altamirano. D., van der Klis, M., Mendez, M., et al. 2008a, ApJ, 685, 436 

—. 2005, Apl, 633, 358 

—. 2008b, Apl, 687, 488 

Belloni, T. & Hasinger, G. 1990, A&A, 227, L33 

Belloni, T., Psaltis, D., & van der Klis, M. 2002, ApJ, 572, 392 

Bradt, H. V., Rothschild, R. E., & Swank, J. H. 1993, A&AS, 97, 355 

Casella, P., Belloni, T„ & Stella, L. 2005, Apl, 629, 403 

Dieters, S. W. & van der Klis, M. 2000, MNRAS, 311, 201 

Fridriksson, J. K. 2011, PhD thesis, Massachusetts Institute of Technology 

Fridriksson, J. K., Homan, J., & Remillard, R. A. 2015, Apl, 809, 52 

Gladstone, J., Done, C., & Gierlinski, M. 2007, MNRAS, 378, 13 

Heil. L. M„ Uttley, P„ & Klein-Wolt, M. 2015, MNRAS, 448, 3348 

Homan, J. 2012, ApJL, 760, L30 

Homan, I., Jonker, P. G., Wijnands, R., van der Klis, M., & van Paradijs, J. 
1999, ApJL, 516, L91 

Homan, J. & van der Klis, M. 2000, ApJ, 539, 847 

Ingram, A. & Done, C. 2011, MNRAS, 415, 2323 

Ingram, A., Done, C., & Fragile, P. C. 2009, MNRAS, 397, L101 

Ingram, A. & van der Klis, M. 2015, MNRAS, 446, 3516 

Jahoda, K., Markwardt, C. B., Radeva, Y., et al. 2006, ApJS, 163, 401 

Jonker, P. G., van der Klis. M., Homan, J., et al. 2000, ApJ, 531, 453 

Jonker, P. G., van der Klis. M., & Wijnands, R. 1999, ApJL, 511, L41 

Lai, D. 1999, ApJ, 524, 1030 

Levine, A. M„ Bradt, H„ Cui, W„ et al. 1996, ApJL, 469, L33 
Miller, J. M. & Homan, J. 2005, ApJL, 618, L107 
Miller, M. C. 1999, ApJ, 520, 256 

Miyamoto, S., Kimura, K., Kitamoto, S., Dotani, T., & Ebisawa, K. 1991, 
ApJ, 383, 784 

Motta, S. E„ Casella, R, Henze, M„ et al. 2015, MNRAS, 447, 2059 
Munoz-Darias, T., Fender, R. R, Motta, S. E., & Belloni, T. M. 2014, 
MNRAS, 443, 3270 

Ponti, G., Munoz-Darias, T., & Fender, R. P. 2014, MNRAS, 444, 1829 

Psaltis, D., Belloni, T., & van der Klis, M. 1999, ApJ, 520, 262 

Remillard, R. A., Lin, D., Cooper, R. L.. & Narayan, R. 2006, ApJ, 646, 407 

Schnittman, J. D., Homan, J., & Miller, J. M. 2006, ApJ, 642, 420 

Shirakawa, A. & Lai, D. 2002, ApJ, 564, 361 

Stella, L. & Vietri, M. 1998, ApJL, 492, L59 

van der Klis, M. 1995, in Proceedings of the NATO Advanced Study 

Institute on the Lives of the Neutron Stars, held in Kernel', Turkey, August 
29-September 12, 1993. Editors, M.A. Alpar, U. Kiziloglu, and J. van 
Paradijs; Publisher. Kluwer Academic. Dordrecht, The Netherlands, 
Boston, Massachusetts, p. 301 

van Straaten, S., van der Klis, M., & Mendez, M. 2003, ApJ, 596, 1155 
Wijnands, R. & van der Klis, M. 1999, ApJ, 514, 939 
Wolff, M„ Ray, P„ Wood, K., & Wijnands, R. 2008a, The Astronomer’s 
Telegram, 1812 

Wolff, M. T., Ray, P. S., & Wood, K. S. 2008b, The Astronomer’s Telegram, 
1736 

Wolff, M. T„ Ray, P. S., Wood, K. S., & Hertz, P. L. 2009, ApJS, 183, 156 




